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A mid-infrared quantum cascade laser with high wall-plug efficiency is analyzed by means of an ensemble 
Monte-Carlo method. Both the carrier transport and the cavity field dynamics are included in the simulation, 
offering a self-consistent approach for analyzing and optimizing the laser operation. It is shown that at low 
temperatures, photon emission and absorption can govern the carrier transport in such devices. Further- 
more we find that photon-induced scattering can strongly affect the kinetic electron distributions within the 
subbands. Our results are validated against available experimental data. 



I. INTRODUCTION 

Quantum cascade lasers (QCLs) are among the most 
promising mid-infrared (MIR) laser sources, offering ap- 
plications in gas sensing, free space communications and 
spectroscopy. Since the first operating device was pre- 
sented in 1994,-*^ QCL designs have been constantly im- 
proved with respect to their efficiency and output power. 
Recently QCLs with wall-plug efficiencies (WPEs) of 
around 50% were reported for the first timei^i^ In such 
structures, light emission and absorption are not only 
relevant with respect to the generated optical power, 
but also strongly affect the carrier transport in the de- 
vices. In fact, for the low temperatures where these 
high efficiencies are reached, the photon-induced pro- 
cesses dominate the other scattering mechanisms. Thus, 
to adequately model the operation of these lasers, the 
optical cavity field has to be considered in the simula- 
tion. While this is routinely done in one-dimensional 
simulations r'"— the cavity field is usually neglected in 
fully three-dimensional (3-D) approaches like the ensem- 
ble Monte-Carlo (EMC)j- non-equilibrium Green's func- 
tions (NEGF)SiiS or 3-D density matri»iiii^ method. 
However, such 3-D simulations, where the in-plane car- 
rier dynamics is explicitly considered, do not only yield 
level occupations, but also the kinetic carrier distribu- 
tions within these levels. Here, we employ the EMC 
method, which has has been intensely used to inves- 
tigate the carrier transport in both MIR^iii^"— and 
terahertaiSr^ QCLs. To include also the optical pro- 
cesses, we have recently extended this approach, allow- 
ing for self-consistent coupled simulations of the car- 
rier transport and the optical cavity field<^ The EMC 
method is a semiclassical approach, i.e., quantum corre- 
lations are neglected in contrast to NEGF-^ or density 
matrijii^ calculations; however, the carrier transport in 
MIR QCLs has been shown to be largely incoherent J^ 



The goal of the present study is to analyze the carrier 
transport and lasing operation in a record-efficiency MIR 
QCL,'^ with a particular focus on the influence of photon- 
induced scattering on the carrier transport. Specifically, 
we show that the inclusion of light emission and absorp- 
tion in the simulation is crucial to obtain a realistic de- 
scription for such devices. Furthermore, our analysis pro- 
vides insight into the carrier dynamics on a microscopic 
level, for example the kinetic electron distributions in the 
upper and lower laser level which are hardly accessible to 
experimental observation. The paper is organized as fol- 
lows: In Sectionini we give a basic overview of our Monte- 
Carlo approach, specifically adapted for InGaAs/InAlAs 
strain-compensated MIR QCLs. In Section Hill we com- 
pare our simulation results to available experimental data 
and demonstrate the strong influence of stimulated emis- 
sion on the carrier transport. The paper is concluded in 
Section HVl 
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II. METHOD 

The EMC method is based on the semiclassical Boltz- 
mann transport equation. ^^ Scattering is self-consistently 
accounted for based on Fermi's golden rule. All the 
relevant mechanisms like electron (e)-longitudinal opti- 
cal (LO) and acoustic phonon, e-interface roughness, e- 
impurity and e-e scattering are routinely considered in 
our simulation too L^^'^^ Moreover, various effects rele- 
vant for MIR QCLs based on the InGaAs/InAlAs mate- 
rial system have been added. We have included random 
alloy scattering,^* with a scattering potential of 0.3 eV re- 
ported for high indium content InGaAs J?^ Furthermore, 
we account for InAs- and GaAs-like phonons, using their 
composition dependent values for the phonon energyj^ 
The scattering rates are weighted by the concentration 
of the individual materials (InAs and GaAs) . The influ- 
ence of the AlAs-like branch is believed to be negligible 
in QCL structures. "^^ Here, the bulk phonon approxima- 
tion is adopted, which was shown to be a valid approach 
for the simulation of such QCL structures J^ 



The (parallel and perpendicular) effective masses 
have been implemented considering strain^ and non- 
parabolicity. Our implementation of non-parabolicity is 
based on the approach developed by Ekenberg.'^^ Non- 
parabolicity parameters were determined from the mate- 
rial bandgap,— using temperature dependent values.-^ 
In the InGaAs material system, the parallel non- 
parabolicity is enhanced by a factor of 1.7 as compared 
to the perpendicular valuci^ The perpendicular effective 
mass affects the subband energies and wavefunctions, as 
considered in our Schrodinger-Poisson solver j^ The par- 
allel effective mass is accounted for by assigning a differ- 
ent value to each subband, affecting the scattering rates 
in the EMC solver. Here we focus on simulations at a lat- 
tice temperature of 40 K where the investigated structure 
operates with a record wall-plug efficiency of above 50% 4 
At such low temperatures, the kinetic electron energies 
are still moderate, whereas for room temperature opera- 
tion, a more complex implementation of nonparabolicity 
might be required, e.g., based on k-p theory J^"— Fur- 
thermore, at low temperatures, the electron leakage into 
indirect valleys, not considered in our simulations, is very 
small Ji.i^ 
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FIG. 1. (Color online) Current-voltage characteristics. The 
EMC simulation results with (X marks) and without (crosses) 
lasing included are compared to available experimental datai^ 
(solid curve). The electric current is governed by stimulated 
photon emission and absorption processes, as has to be ex- 
pected for a WPE as high as 50%. 



III. RESULTS 



The interface roughness is typically described by a 
mean height A and a correlation length A. In contrast to 
the well-known bulk material parameters, this quantity is 
hardly accessible to experimental measurement and de- 
pends critically on the growth conditions. Thus, there is 
an uncertainty regarding the values of A and A.^iS2 How- 
ever, experimental data indicate that AA w 1 nm^ for the 
InGaAs/InAlAs structures y^^^ reducing the uncertainty 
to a single parameter value. We choose A = 0.06 nm, 
which yields the best agreement with the experimental 
results. This value is somewhat lower than previously 
used values for strain-free lattice-matched structures i^i^ 
However, we note that vertical correlations, which are 
not included in our simulations, can reduce the effect 
of interface roughness for strained (e.g., strain-balanced) 
quantum cascade lasers,- - as considered here. 



Lasing is implemented based on a recently published 
approach, treating the photon dynamics in terms of clas- 
sical intensity evolution equations and accounting for 
photon-induced scattering in the EMC solver i2M^ In this 
way we can self-consistently describe the coupled carrier- 
light dynamics due to absorption as well as stimulated 
and spontaneous emission. For the investigated design 
operating at 5 /im,'^ the mirror loss, which amounts to 
6.4 cm~^ for a 2 mm long structure, dominates the waveg- 
uide loss, which is about 0.5 cm^^ for such cavities 4^ The 
confinement factor is chosen to be 0.8 as found for a simi- 
lar design J^ For our simulation, we use 1200 longitudinal 
modes in the frequency range between 50 and 80 THz, 
corresponding to a Fabry-Perot mode spacing of 25 GHz. 



Results are presented for a recently fabricated high ef- 
ficiency QCL operating at 5/^m.'^ The simulations were 
performed at a lattice temperature of 40 K, where the 
record WPE of 53% was observed. 

In Fig. [U we compare the current-voltage charac- 
teristics provided by EMC simulations to experiment. 
The simulations were performed at biases ranging from 
115kV/cm to 145kV/cm. For comparison to experi- 
ment, these were converted to the voltage points in Fig. [I] 
by considering 80 stages with a thickness of 22.1 nm 
each,-^ Above threshold, good agreement is found if lasing 
is included, while the current due to non-radiative pro- 
cesses (EMC without lasing) is lower by a factor of almost 
3 than the experimentally measured current. This shows 
that stimulated processes become more and more impor- 
tant for a correct description of the carrier transport as 
the WPE of QCLs is improved. On the other hand, the 
spontaneous photon emission rates in our simulation are 
far too low to affect the carrier transport, which is in 
agreement with theoretical considerations.™ The onset of 
the negative differential resistance (NDR) regime agrees 
well, occuring at 25.1 V for the EMC with lasing included 
and 25.6 V in the experiment. For low fields where the 
energy levels are not aligned, the simulation underesti- 
mates the experimentally observed current. Here, the 
scattering- induced transport is not efficient, and the re- 
maining current can likely be attributed to coherent low- 
field transport which is not included in the EMC simula- 
tion!^ For design optimization with respect to the WPE, 
the parasitic channels should be suppressed and the stim- 
ulated emission into the lasing modes maximized. Such 
a task can only be performed with an approach taking 
into account the optical cavity field. 

In Fig. [Ifa) and (b) we compare the simulated and ex- 
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FIG. 2. (Color online) (a) Current-optical power and (b) 
current- WPE characteristics. The EMC simulation results 
with lasing included (X marks) are compared to available ex- 
perimental data- (solid curve). 



perimental current-output power and current- WPE char- 
acteristics. In the EMC simulation, the bias dependent 
WPE r/wPE is computed as ?7wpe = Po-pt/Pc\- Here, 
^opt is the simulated optical power emitted through both 
facets as in the experiment,"^ and the electric power Pel is 
the product of the applied voltage and the simulated elec- 
tric current. The simulated and experimental current- 
output power characteristics in Fig. ^a) show excellent 
qualitative and quantitative agreement. The maximum 
emitted optical power is about 10 W, which is in both 
cases obtained around the onset of NDR, where the cur- 
rent reaches its maximum value of 0.8 A. For higher bi- 
ases, i.e., in the NDR regime, the simulated optical power 
and electric current get reduced again. The simulated 
threshold current is lower than the experimental value, 
for the reasons discussed in the previous paragraph. Also 
the simulated and experimental current- WPE character- 
istics shown in Fig. Hl^b) agree well. Particularly, the 
maximum simulated WPE of 49% below the onset of 
NDR compares very well to the experimental value of 
53%. The simulated high WPE value of 53.5% around 
the onset of NDR is not observed in the experiment, 
which we attribute to the fact that the operation in the 
NDR region is unstable due to domain formation4^i^ 

Full fc-space three-dimensional simulation approaches 
like EMC can yield information on the microscopic level, 
which is hardly accessible to experimental observation. 
In the following, we investigate the intrasubband ki- 
netic carrier distributions. These can be characterized 
by corresponding electron temperatures only in the case 
of quasi-thermal equilibrium within the subbands, cor- 
responding to a Maxwellian distribution for low doping. 
The kinetic electron distribution in the upper and lower 
laser level is shown in Fig.[31[a). The bias is 25.1 V, where 
the simulated current and output power reach their max- 
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FIG. 3. (Color online) (a) Simulated kinetic electron dis- 
tribution / {E) in the upper and lower laser level with and 
without lasing included, (b) Energy resolved electron density 
without lasing included; (c) energy resolved electron density 
with lasing included. The upper laser levels are marked by 
solid rectangles, and the lower laser level is marked by dashed 
rectangles. 



imum values. We note that for this bias, optical transi- 
tions from two upper levels contribute significantly to 
lasing. For simplicity, we restrict our discussion to one of 
these states, since the kinetic electron distribution func- 
tion is found to be similar for the other level. The carrier 
distributions in the laser levels change significantly by 
taking lasing into account (solid curves), as compared to 
the case without lasing (dashed curves). The lasing ac- 
tion leads to a depletion of the upper laser levels and a 
filling of the lower laser level, corresponding to the ef- 
fect of gain saturation. In the inset of Fig. [2l^a), the 
electron distribution in the lower laser level without las- 
ing is shown on a logarithmic scale, i.e., a Maxwellian 
distribution would produce a straight line. The distri- 
bution is highly non-Maxwellian with an additional peak 
at around 250 meV, corresponding to the energy spacing 
between upper and lower laser level. This bump stems 
from nonradiative transitions from the upper laser level, 
mainly LO phonon scattering as the dominant nonradia- 



tive mechanism. The energetic extension of the bump is 
partly due to the kinetic electron distribution in the up- 
per laser level and the finite phonon energies of 29.5 meV 
and 32.2 meV for the two LO branches considered here. 
E-e scattering, which is the predominant intrasubband 
scattering mechanism, is unable to thermalize the intra- 
subband carrier distribution strongly enough to suppress 
the bump;^!^ The inclusion of lasing action leads not 
only to a filling of the lower laser level, but also to a 
more thermalized kinetic electron distribution, while the 
bump at around 250 meV still persists. A least square fit 
produces an electron temperature of Te = 314 K (upper 
laser level) and Te = 344 K (lower laser level) with lasing 
included. This is consistent with the observation that in 
strain compensated structures the electronic temperature 
is clearly above the lattice temperature i^i^ 

In Fig.[3Jb) and (c), the energy resolved electron den- 
sity (normalized to its maximum value) is shown without 
and with lasing included, again for a bias of 25.1V. For 
the two quantum wells located between 6 nm and 18 nm, 
the upper laser level is omitted, so that the high-energy 
peak discussed in the previous paragraph (see inset of 
Fig. |3Ka)), here located at around 1.9 eV, is clearly vis- 
ible. By comparing Fig. |3Kb) and (c), we can observe 
the changes in the energy resolved electron density for 
the laser levels without and with lasing included. The 
high-energy tails of the kinetic electron distributions re- 
main basically unaffected. In particular, the additional 
high-energy peak for the lower laser level appears also if 
lasing is accounted for. 

We have successfully tested our approach for a further 
high efficiency QCL design^ using the same material pa- 
rameters, again finding good agreement with experiment. 



IV. CONCLUSION 

Based on a self-consistent EMC carrier transport simu- 
lation including the optical cavity field, we have analyzed 
the effect of photon emission and absorption on the car- 
rier transport in a high WPE quantum cascade laser. 
In the regime where efficient lasing is obtained, we find 
that the inclusion of photon-induced scattering is crucial 
for the correct calculation of the device current. Fur- 
thermore, a comparison to experimental data yields very 
good agreement for the optical output power and WPE. 
An analysis of those quantities, as also needed for design 
optimization, is only possible with an approach which in- 
cludes both the carrier transport and the optical cavity 
field. The EMC method also enables us to investigate 
microscopic quantities such as the intrasubband kinetic 
carrier distributions, hardly accessible to experimental 
observation. The upper and lower laser level carrier dis- 
tributions are strongly affected by the lasing action and 
approach each other, corresponding to gain saturation. 
We observe strong deviations from an equilibrium distri- 
bution especially for the lower laser level, where a high- 
energy peak in the electron distribution is found, caused 



by parasitic transitions from the upper laser level. Our 
results show that the chosen approach is well suited to 
model high efficiency MIR QCLs on a qualitative and 
quantitative level, and to analyze the laser operation on 
a microscopic scale which is hardly accessible to experi- 
mental observations. 
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